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The question whether mean held effects exist in heavy-ion collisions at AGS 
energies is studied in the framework of A Relativistic Transport (ART) model. 
It is found that in central collisions of Au-|-Au at Pbeam/A =11.6 GeV/c a simple, 
Skyrme-type nuclear mean field satisfying the causality requirement reduces the 
maximum baryon and energy densities reached in the cascade model by about 
30% and 40%, respectively. The mean field affects the inclusive, single particle 
observables of various hadrons by at most 20%. We show, however, that the 
mean field causes a factor of 2.5 increase in the strength of the baryon transverse 
collective flow. 
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The purpose of relativistic heavy-ion collisions at AGS energies is to study the properties 
of hot and dense hadronic matter and the possible phase transition to the quark-gluon 
plasma. A vast body of data have been collected and analysed during the past years [|^. 
Comparisons of these data with theoretical models, such as RQMD [Q, ARC [Q and QGSM 
1^, have revealed much interesting physics. In particular, a picture of nearly complete 
stopping of baryons in central heavy-ion collisions at AGS energies has emerged from these 
studies. It has been shown that baryon and energy densities up to 10 times that of nuclei in 
their ground state have been reached in these collisions. This has led to the suggestion that 
the energy and baryon densities achieved during central heavy-ion collisions at AGS energies 
may be high enough to form the quark-gluon plasma 0. However, the above suggestion 
and many speculations critically rely on the inferred maximum energy and baryon densities 
reached in the reaction. 

In this Letter, we report on the main aspects of a relativistic transport (ART) model for 
heavy-ion collisions at AGS energies, and study, in particular, effects of the nuclear mean held 
which has been ignored in ARC. Our motivations of studying the nuclear mean held ehects 
are mainly the following. Firstly, we believe that the mean-held potential is not negligible in 
heavy-ion collisions at AGS energies. Although the forward scattering amplitudes of hadron- 
hadron collisions in the high energy limit have been found approximately purely imaginary 
@], the AGS energies may not be high enough for the real part of the scattering amplitude 
to completely vanish. Of course, the form and strength of the corresponding mean held 
in the hot and dense hadronic matter is highly uncertain and has been a subject of much 
discussions. Secondly, although the kinetic energy is much higher than the potential energy 
in the early stage of the reaction, particles are gradually slowed down and the mean held 
plays an increasingly important role as the reaction goes on. In particular, the repulsive 
mean held in the high density region tends to keep particles from coming too close to 
each other and therefore reduce the maximum energy and baryon densities reached in the 
reaction should there be no mean held. Moreover, in the expansion phase of the reaction 
mean-held ehects are expected to be even stronger. It is therefore necessary to study how 
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the compression and expansion are affected by the nuclear mean field. Finally, what are the 
experimentally observable consequences of the nuclear mean field ? 

Our relativistic transport model is developed from the well known Boltzmann-Uhling- 
Uehlingbeck (BUU) model (e.g. []^,^) for intermediate energy heavy-ion collisions. Here we 
report on the extension of the model to AGS energies with a reasonable amount of details. We 
have included in the model the following baryons: N, A(1232), A*(1440), A*(1535), A, S, 
and mesons: tt, p, cu, p, K, with their explicit isospin degrees of freedom. Both elastic 
and inelastic collisions among most of these particles are simulated as best as we can using 
as much input from the experimental hadron-hadron data as possible |^. Most of the 
inelastic hadron-hadron collisions are modeled through the formation of resonances. The 
finite lifetimes of these resonances effectively take into account the formation time of newly 
produced secondaries. More specifically, we have included in the model 


NN^NA, AA*(1440), AA*(1535), 

(1) 

NN ^ A A, AA*(1440), 

(2) 

NN NNp, NNu, AAvr, 

(3) 

NN AAp, 

(4) 

NA ^ AA*(1440), AA*(1535), 

(5) 

AA AA*(1440), AA*(1535), 

(6) 

AA*(1440) ^ AA*(1535), 

(7) 

and those producing kaons 

NN ^ na{j:)k, aa(s)a:. 

(8) 

NR aa(s)a:, aa(s)a:. 

(9) 

rr^aa(s)a:, aa(s)a:. 

(10) 


where R is A, A*(1440) or A*(1535). Cross sections parameterized by VerWest et ah M 
are used for the single A and A* (1440) production in processes (1). The cross section of 
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A^*(1535) production is estimated from that of 7] production |TT[. Cross sections for the 
double resonances production in processes (2) are taken to be the same for all possible 
isospin channels, this is supported by recent calculations based on the one-boson-exchange 


model 1]^. Numerically, it is estimated by subtracting from the inclusive 27r production 
cross section the contribution of NN NNp and the 27r decay of the iV*(1440) in 
the NN —> iViV*(1440) process. The cross sections for p and oo production in channels 
(3) are taken directly from the experimental data |^. The cross section for the process 
NN —>■ AAvr is taken as the difference between the inclusive dvr and the u production 
cross section. We attribute the difference between the experimental total nucleon-nucleon 
inelastic cross section and the sum of cross sections for channels (1) to (3) as well as the 
kaon production cross sections of channels (8) to the process NN AAp. This is done 
to ensure the total inelasticity of nucleon-nucleon collisions as only a limited, though large, 
number of reaction channels have been incorporated so far. The bias introduced by this 
approximation towards the cross section of the quasi-dvr production process NN —>■ AAp is 
very small. Cross sections for channels (5) to (7), (9) and (10) are taken to be the same as 
in nucleon-nucleon collisions having the same center of mass energy and total charge. Cross 
sections for the inverse processes are calculated by using the detailed balance. Masses of 
baryon and meson resonances are generated according to the single or joint Breit-Wigner 
distributions with momentum dependent widths using the rejection method. 

One can also separate meson-baryon collisions into elastic and inelastic parts. We 
model elastic collisions through both the formation of baryon resonances, i.e., vriV <->■ 
A(A*(1440), A*(1535)) and rjN <->■ W(1535), as well as direct processes, i.e., 7i{p) + 
N{A, N*) —>• 7r(p) -|- N{A, N*), and K + N ^ K + N. The formation of the three baryon 
resonances accounts almost completely the vr -|- A elastic cross sections at low energies. At 
higher energies, for 7r“ -|-p at ^/s > 1.7GeV, for example, the elastic cross section is about 
7 mb and is mainly due to the formation of higher resonances which are not included in 
the present model. We therefore attribute the difference between the experimental elastic 
cross section and the contribution from the three baryon resonances to the direct process 
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7T + N —> vr + A^. For experimentally unknown cross sections, such as, + 7t + A{N*) and 
p + N, we calculate them using the resonance model. Neglecting the interference between 
resonances, one has 




(2J^ + 1) 


TUM + B) 


k'^ ^ {2 Sm + l)( 2 S's + 1) (v^ - + 0 . 25 r|(toto/) 


( 11 ) 


The factor 1.3 is obtained by htting to the high energy part of the 7r’''+p data. For calculating 
the cross sections of direct channels, R runs over all baryon resonances with masses up to 2 
GeV. 

The 71 + N inelastic collision mainly goes through the production of pions and kaons. 
Similar to the treatment of baryon-baryon collisions, we model the inelastic vr + iV 
collisions through the production of resonances, namely, vr + A + 7r{p,uj) and 

'k{p,uj) + N{A,N*) —>■ A(E) + K. All experimental cross sections for tt + A collisions 
with hnal states having two and three pions are attributed to the production of Avr and Ap, 
respectively. The difference between the experimental total vr + A inelastic cross section 
and cross sections for the production of two and three pions as well as kaons are attributed 
to the production of Ao;. Cross sections for the production of kaons in 7r(p, a;) + A(A, A*) 
collisions are taken to be the same as that in the vr + A collision. 

We model pion-pion elastic collisions through the formation of the p meson and the direct 
process. The latter takes into account the case when the quantum numbers of colliding pions 
forbid the formation of the p meson. These cross sections are taken mainly from ref. |]T5 


The inelastic collisions among mesons are modeled through the production of KK. The 
cross section for this process is highly uncertain, here we use the cross section calculated 
from the A*-exchange model of ref. H- 


The “leading” particle behaviour in energetic collisions is ensured by requiring the out¬ 
going baryons to have the same or similar quantum numbers as the incident ones so that 
their longitudinal directions are retained when performing the momentum transformation 
from the baryon-baryon c.m. frame to the nucleus-nucleus c.m. frame. The longitudinal 
and transverse momenta of baryons in the hnal state is assigned according to the following 
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distribution under the constraint of total energy and momentum conservations 
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d'^N 9 9 

—^ cx (1. + O.Sx* - 0.9a;* + 0.5e-^°P*), (12) 

d^p* 

where x* is the scaled longitudinal momentum in the c.m. frame, i.e., x* = 2p*^j\fs. The 
distribution has the properties of the naive scaling, pt and p* factorization and the soft, 
energy independent transverse momentum distribution. With this distribution and the cross 
sections discussed above, we can well reproduce the longitudinal and transverse momentum 
distributions of protons and pions in pp collisions at pheam =12 GeV/c and below. The test 
against pp data together with many other details of the model will be published separately 

0 . 

In the framework of the ART model outlined above, we shall concentrate on studying the 
effects of nuclear mean held. Without much reliable knowledge about the nuclear equation 
of state in hot and dense medium, we use here the simple, Skyrme-type parameterization 
widely used at Bevalac energies and below. 


f/(p) = _358.1(-^) + 304.8(—(13) 
Po Po 

This is the so-called soft equation of state which satishes the causality requirement upto 
about 7po- The standard test-particle method is used to calculate the global baryon density 
Pg and energy density Cg in the nucleus-nucleus c.m. frame on a lattice with the cell size of 
1 fm^. The local baryon density and energy density in each cell are obtained by p; = 
and Cl = e^/y, where y is the Lorentz factor in each cell. 

Let us hrst study effects of the mean hied on the creation of high baryon and energy 
densities. We show in Fig. 1 the evolution of the local baryon and energy densities in 
the central cell during the reaction of Au-|-Au at Pbeam/A=11.6 GeV/c and b=0 for three 
diherent cases. For the cascade case, it is seen that the maximum baryon density of about 
9po and the maximum energy density of about 3.6 GeV/fm^ are reached at about 4 fm/c. 
The high energy density matter lasts for about 5 fm/c. A more detailed examination of the 
density contour plots shows that the high energy density region with ei > 2.0 GeV/fm^ has 


6 







the maximum volume of about 200 fm^ at about 4 fm/c. It is worth mentioning that the 
main features of the cascade model calculations are very similar to that of ARC. 

The currently estimated critical baryon and energy densities for forming the quark gluon 
plasma is about 5po and 2.5 GeV/fm^, respectively p. According to the above cascade 
model prediction, it is highly possible to form the quark-gluon plasma at AGS energies. 
However, the lifetime and volume of the high density region are signihcantly reduced by the 
repulsive mean held. Noticeably, the maximum baryon and energy densities are reduced to 
about 7po and 2.6 GeV/fm^, respectively, by using the soft nuclear equation of state. With 
a stiff equation of state corresponding to the compressibility of K = 380 MeV the reduction 
is even larger although the mean held has almost no ehect in the early stage of the reaction 
when the kinetic energy is much higher than the potential energy. Since the stih equation 
of state violates causality already at about 3po, we will only use the soft equation of state 
in the following. The reduction of the maximum baryon and energy densities due to the 
mean held may be large enough to ahect signihcantly the kind of physical processes that 
can happen during the reaction. 

The validity of any hadronic model for heavy ion collisions at AGS energies is naturally 
limited by the possible phase transition to the quark-gluon plasma. With this precaution 
in mind, one can use predictions of hadronic models as a baseline in searching for new 
phenomena. The formation of the quark-gluon plasma or a mixed phase of hadrons and 
quark-gluon plasma is expected to reduce the pressure of the system and leads to a softened 
nuclear equation of state. It is therefore interesting to search for experimental observables 
that are sensitive to the nuclear equation of state. A careful study on the transverse mo¬ 
mentum and rapidity distributions of protons, pions and kaons for the reaction of Au-|-Au 
at Pfoeam/A =11.6 GcV/c has been done by comparing calculations with and without the 
mean held. Much to our surprise, both calculations can well reproduce the inclusive, single 
particle data within current systematic and statistical error bars of the data. The largest 
difference of about 20% is found for the proton transverse momentum distribution with a 
slight shift to higher transverse momenta in the case using the soft equation of state. This 
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difference is, however, compatible with the current systematic error bars of the data. These 
results will be published elsewhere |n . 


Hinted by the hndings at Bevelac energies that the collective variables or correlation 
functions, unlike the single particle observables, are very sensitive to the nuclear equation of 
state, we now turn to the analysis of the baryon transverse collective flow [I^. Here we use 
the standard flow analysis as in ref. 0. It was shown recently in ref. [Q that an improved 


flow analysis with the explicit conservation of reaction plane and angular momentum in 
individual hadron-hadron collisions results in an increase of the transverse momentum by 
about 8% to 23% at Bevalac energies. This small enhancement will not affect our discussions 
and conclusions in the following. 

We show in Fig. 2 the average transverse momentum of nucleons in the reaction plane 
as a function of rapidity for the reaction of Au-|-Au at Pfeeam/A =11.6 GeV/c and impact 
parameters of 2, 6 and 10 fm. Signihcant differences exist between calculations with and 
without the mean held for the reaction at all of the three impact parameters. In partic¬ 
ular, the how parameter dehned as the slope of the transverse momentum distribution at 
midrapidity is about a factor of 2.5 larger in the case with the mean held. The strength of 
the so-called “bounce-oh” ehect at target or projectile rapidities is also much stronger in 
calculations with the mean held. It is also seen that the collective how is the strongest in 
the midcentral collisions. 

To see why the collective how is a sensitive probe of the dynamics in the high density 
region, we show in Fig. 3 the how parameter as a function of time for the Au-|-Au reaction 
at impact parameters of 2 fm and 6 fm. It is clear that the how is mainly generated in the 
high density region and saturates in the expansion phase. The ratio of hnal how parameters 
in the calculations with and without the mean held is about 2.5 in the reactions at both 
impact parameters. This diherence is large enough for a clear distinction between models 
when the experimental data is available. For the reaction at an impact parameter of 6 fm, 
the how parameter undergoes a reduction before saturation. This is due to the rehection 
of hot baryons from the cold spectator nucleons. It is interesting to note that the decrease 
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of pressure in the high density region due to the possible phase transition to guark-gluon 
plasma is expected to reduce the strength of collective flow. An experimentally measured, 
signihcantly smaller flow parameter than the cascade prediction will be a strong indication 
of the quark gluon plasma formation at AGS energies. 

In summary, we have developed a new relativistic transport model for heavy-ion collisions 
at AGS energies. Within the framework of this model, we have found that the mean field 
can signihcantly affect the maximum energy and baryon densities reached in the reaction. 
The transverse collective how analysis is suggested as a useful tool to disentangle theoretical 
models and a possible indicator of the formation of quark-gluon plasma at AGS energies. 
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Koch, J. Randrup and Gy. Wolf for their constructive suggestions and criticisms. This 
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FIGURE CAPTIONS 

Fig. 1 Evolution of the central local baryon density and energy density in the reaction 
of Au-|-Au at Pbeam/A. =11.6 GeV/c and b=0. 

Fig. 2 Baryon average transverse velocity in the reaction plane as a function of rapidity 
for Au-|-Au reactions at Pbeam/A=11.6 GeV/c and three different impact parameters. 

Fig. 3 Time evolution of the flow parameter in the reaction of Au-|-Au at P6eam/A=11.6 
GeV/c and the impact parameter of 2 and 6 fm. 
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